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ABSTRACT 

We present evidence for the existence of significant substructure in the cold front cluster Abell 3667 
based on multi-object spectroscopy taken with the 3.9m Anglo Australian Telescope. This paper is the 
second in a series analyzing the relationship between cold fronts observed in Chandra X-ray images 
and merger activity observed at optical wavelengths. We have obtained 910 galaxy redshifts in the 
field of Abell 3667 out to 3.5 Mpc, of which 550 are confirmed cluster members, more than doubling the 
number of spectroscopically confirmed members previously available and probing some 3 mag down the 
luminosity function. From this sample, we derive a cluster redshift of z = 0.0553 ±0.0002 and velocity 
dispersion of 1056 ± 38kms _1 and use a number of statistical tests to search for substructure. We 
find significant evidence for substructure in the spatial distribution of member galaxies and also in the 
localized velocity distributions and, in spite of this evidence, find the global velocity distribution does 
not deviate significantly from a Gaussian. Using combined spatial and velocity information, we find 
the cluster can be separated into two major structures, with roughly equal velocity dispersions, but 
offset in peculiar velocity from each other by ~ 500 km s -1 , and a number of minor substructures. We 
propose two scenarios which explain the radio and X-ray observations. Our data show the cold front 
is directly related to cluster merger activity, and also highlights the extent of optical data required to 
unambiguously detect the presence of substructure. 
Subject headings: galaxies: clusters: individual (Abell 3667) 



1. INTRODUCTION 

At the pinnacle of the structure formation hierarchy, 
major mergers of clusters involve the largest, most mas- 
sive virialized objects at the present epoch. In a major 
merger, up to 10 64 erg of gravitational potential energy 
is dispersed among the cluster constituents, with around 
15% going into the gas through shock heating, dissipa- 
tion of turbulence and adiabatic compression of the in- 
tracluster medium (ICM). The majority of the energy 
goes into inducing peculiar motions of the galaxies and 
dark matter, and a small fraction into cosmic ray accel- 
eration and magnetic field amplification. Despite some 
evidence pointing to rapid evolution of gal axies caused 
by a major merger (jCaldwell fc Rosdll997h . the impli- 
cations of mergers for galaxy evolution are not well un- 
derstood, partly because of the difficulty in identifying 
major mergers, and also in identifying galaxies belong- 
ing to substructure. To disentangle the physics behind 
each of these merger driven processes in the ICM, and 
also to robustly identify substructure within clusters, a 
multi-wavelength approach is required. 

To this end, we have selected a sample of clusters ex- 
hibiting robust examp les of cold fronts from the Chandra 
archive (jOwers 1120081 ). and this is the second in a series 
of papers exploring the sample by emplo ying a multi- 
wave length analysis. In the first paper (jOwers et al.l 
2009, hereafter Paper I), we present a combined X- 
ray and optical analysis of Abell 1201, finding multiple 
substructures through combination of redshift and spa- 
tial information for 321 cluster members and confirming 
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that the cold fronts observed are a direct result of on- 
going merger activity. Here, we turn our attention to 
Abell 3667 which is an integral component of this ongo- 
ing study, gi ven it was one of the firs t cold fronts detected 
by Chandra (|Vikhlinin et al.ll2001af ) and helped to define 
the cold front phenomenon. Th e largest spectroscopic 
study to date was carried out by iJohnston-Hollitt et alJ 
(2008) and based on 231 spectroscopically confirmed 
cluster members. They find no evidence for dynami- 
cal activity in the velocity distribution, nor any evidence 
of localized deviations from the global velocity structure 
and conclude that in order to rule out the existence of 
substructure, there is a need for a larger sample of cluster 
member redshifts. The current study uses more extensive 
optical data for Abell 3667, with the aim of unambigu- 
ously detecting dynamical substructure in Abell 3667 and 
relating it to the cold front. 

At X-ray wavelengths Abell 3 667 is luminous (Lx (0-4— 
2.4 keV) = 5.1 x lO^ergs" 1 : labeling et al.l ll996h. hot 
(kT ~ 6 keV; lKnopp et al.lll996l ) and displays several fea- 
tures indicative of a major merger. The ROSAT X-ray 
observations of iKnopp et alJ (|1996l ) suggested the clus- 
ter was not in dynamical equilibrium, with an overall 
elliptical appearance, central structure and significant X- 
ray emission associated with the second D galaxy to the 
north-west. Perhaps the most interesting feature in the 
X-ray images is the surface brightness discontinuity to 
the south-east of th e cluster center initially interpreted 
as a merger shock bv lMarkevitch et al.l (|1999h . However, 
the increase in both spectral and spatial resolution and 
sensitivity with the Chandra X-ray observatory provided 
the necessary data to conclude that the discontinuity is 
a cold front - a contact discontinuity at the interface of 
dense, cold gas cloud movin g at near sonic speed th rough 
the hotter, less dense ICM ijVikhlinin et alil2001ah . The 
striking sharpness of the front suggests transport pro- 
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cesses across the front must be suppressed, and the sta- 
bility of the front has led to interesting insights into the 
growth of hydrodynamic instabil ities due to the motion 
of the gas cloud within the I CM (jVikhlinin et alj|2001bt 
IChurazov fc Inogamovll2"004h . 

Abell 3667 exhibits more interesting signatures of an 
ongoing merger at radio wavelengths, and distinguishes 
itself as one of only a few clusters with two diffuse ra- 
dio relics positioned symmetrically on opposing sides of 
the cluster periphery with their major a xis perpendicular 
to th e direction of cluster elongation (jRottgering et al.1 
1997). Radio relics are steep sp ectrum, filamentary, po - 
larized, synchrotron sources and lRottgering et al.1 (fl997l ) 
hypothesized that turbulence in the cluster outskirts as- 
sociated with a cluster merger was responsible for re- 
accelerating electrons to the relativistic speeds needed 
to produce the observ ed synchrotron emission. However, 
iRoettiger et al.l (|1999|) were able to produce the observed 
properties of the radio relics by simulating a cluster 
merger which produced shock fronts that inject the en- 
ergy required to accelerate electrons to ultra-relativistic 
speeds required for the synchrotron radiation (also see 
iHoeft fc Brug genl I2007L and references therein). More 
indirect evidence for a merger comes in the form of a 
head-tail galaxy loca ted ~ 635 kpc to the north-west 
of the cen tral galaxy (|Rottgering et alJll997h . which ac- 
cording to Blit on et al.l ( 19981 ) are preferentially found in 
dynamically active clusters. 

In the optical, Abell 3667 manifests itself as a rich- 
ness class 2, Bautz -Morgan I-II intermedia te type cluster 
" Abell et al.lll989h at redshift z = 0.055. iGirardi et all 
1998) measure a virial mass of 1.2 x 10 15 /i _1 M Q and 
velocity dispersion 971 kms for 154 members w ithin 
2.2 fr" 1 Mpc Earli er work by ISodre et all (| 1992ft and 
iProust et al.l (|1988l ) measure higher dispersions (1100- 
1400 km s -1 ) but with smaller samples. Despite the 
clear evidence of merger activity at X-ray and radio 
wavelengths, the existence of substructure detected at 
optical wavelengths remains ambiguous. Multimodal- 
ity has been observ ed in projected gala xy density maps 
(jSodre et all Il992t IProust et all 11988ft . although this 
was not detected at high significance in the density 
maps of spectroscopica l ly con firmed cluster members in 
I Johnston- Hollitt et al.l (|2008l ). B ased on the DEDICA 
substructure detection algorithm, iRamella et al.1 (|2007ft 
detect a number of substructures i n Abell 3667 using spa- 
tial plus color information, whilst [Joffre et afl (|2000ft use 
weak lensing to produce a projected mass map which is 
strongly peaked around the central D galaxy, with three 
other peaks at lower significance - one to the north, one 
coincident with the second bright D galaxy and one to 
the south-east of the cluster center. 

In this paper, we present multi-object spectra (MOS) 
obtained with AAOmega at the 3.9m Anglo Australian 
Telescope (AAT), more than doubling the number of 
spectroscopically confirmed members of Abell 3667. In 
$2] we present the photometric catalog, selection crite- 
ria for MOS observations and data reduction. In Sj3jwe 
present the method for cluster membership determina- 
tion. We outline the methods used for detection of sub- 
structure in Sj4]and in Sj5]we discuss and summarize our 
analysis. Throughout the paper, we assume a standard 
ACDM cosmology where Hq = 70 kms -1 , Q, m = 0.3 and 
Oa = 0.7. For this cosmology and at the redshift of the 



cluster 1" = 1.075 kpc. 

2. OPTICAL SAMPLE SELECTION AND OBSERVATIONS 

In this section, we present the method for selection of 
photometric candidates for spectroscopic follow up, de- 
tails of the observations conducted at the 3.9m Anglo- 
Australian Telescope (AAT) and brief descriptions of 
data reduction, redshift determination, redshift accuracy 
and spectroscopic completeness. 

2.1. Parent Photometric Catalog 

The parent photometric catalog was obtained from the 
Supercosmos Sky Survey (SSS) server 4 using the object 
catalog extraction function. Objects classified as galaxies 
(SSS class=l) within a radius R = 53' (3.4 Mpc) of the 
central dominant cluster galaxy at RA=20 h 12 m 27.38 s , 
DEC=-56° 49' 35.7" were included. We further filtered 
the catalog such that for R < 2 Mpc, only objects with 
magnitudes rp < 19.0 were included, whereas for 2 < 
R < 3.4 Mpc a brighter magnitude limit of rp < 18.0 was 
adopted. These selection criteria limited the catalog to a 
manageable size, given the available telescope time, while 
retaining the objects most likely to be cluster members. 
This resulted in a total of 2,163 galaxies being included 
in our parent catalog. 

To maximize the efficiency of our spectroscopic obser- 
vations in terms of targeting cluster members, we ranked 
each galaxy in the parent catalog on the basis of its po- 
sition in the color-magnitude (CM) plane (see Figure [I}, 
as well as its projected distance from the cluster center. 
The position of galaxies relative to Abell 3667's red color- 
magnitude sequence was used as a means to weight our 
target selection, with the red sequence being identified 
using 160 spectroscopically confirmed members obtained 
from the NED catalog. These, and further spectroscop- 
ically confirmed members as determined in § [31 are rep- 
resented by the green points in Figure [TJ and at magni- 
tudes fainter than rp ~ 16.5 are seen to delineate a red 
sequence which exhibits the usual negative slope in the 
CM diagram. At rp < 16.5, this red sequence persists, 
but exhibits a positive slope, which we can only assume 
is due to some systematic problem in the photometry 
at these brighter magnitudes. This behavior, however, 
is not problematic for pre-selecting targets, since all we 
want to do is use the red sequence to define a likely upper 
red envelope, below which most cluster members reside 
(galaxies above are almost certainly more distant non- 
members). The solid straight line shown in Figure [1] (for 
which bj — rp — 0.168 x rp— 1.771) represents the red en- 
velope that was adopted. Its position and slope were set 
to best include all galaxies on and bluewards of the red 
sequence, allowing for the systematic change in its slope 
and the increasing photometric scatter at fainter magni- 
tudes. Galaxies meeting these criteria at R < 500 kpc 
from the cluster center were given the highest rank. For 
more remote galaxies, the rank was decreased for each 
additional 500 kpc step in distance from the cluster cen- 
ter. Galaxies lying lying above the red envelope were 
all given lower rank than those below, as well as being 
ranked in a similar manner by their cluster-centric dis- 
tances. 

4 see: http: /7www-wfau.roe.ac.uk/sss/| 
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Fig. 1. — Color-magnitude (CM) diagram for objects in the 
field of Abell 3667. The green points represent spectroscopically 
confirmed cluster members (see f[3jl, the dark blue stars are fore- 
ground stars, the light blue stars are foreground galaxies, and the 
red squares are background galaxies. The black points are objects 
without spectra. The solid line shows the red envelope we adopted 
for separating objects which are redward of the cluster's red CM 
sequence and hence are likely non-members, from those that are on 
or blueward of the red sequence and hence have a high likelihood 
of being members. 

For completeness, we also plot in Figure Q] objects 
which have been shown in § to be foreground stars 
(shown as the dark blue stars), foreground galaxies (light 
blue stars), and background galaxies (red squares). The 
small black points indicate objects in the field without 
spectra. 

2.2. AAOmega Observations 

Substantial new spectroscopy of Abell 3667 was ob- 
tained with the Two Degree Field(2dF)/AA0mega 
multi-object spectrograph facility on the 3.9 m Anglo- 
Australian Telescope (AAT) during the nights of 2007 
July 14-18. AAOmega is a bench-mounted dual-beam 
spectrograph which is fed by 400 fibers which are robot- 
ically placed wi thin the AAT's two degree field of view 
at prime- focus (ISaunders et al.ll200l ISmith et all 120041 : 
ISharp et alJ 120061 ). All observations were taken using 
the medium resolution (R ~ 1300) 580V (blue arm) 
and 385R (red arm) gratings. In combination with the 
2 arcsec diameter fibers, they deliver a spectral resolution 
at the detector (a 2kx4k 15 fim pixel E2V CCD in each 
arm) of 3.6A in the blue and 5.5A in the red. The wave- 
length range covered by the two arms is 3700 - 8800A. 

The large number of galaxy targets and their highly 
clustered distribution on the sky, together with 30 arcsec 
restriction on the minimum separation between fibers, 
meant that six different plate configurations were re- 
quired to obtain a high level of spectroscopic complete- 
ness and adequate coverage in the central regions of 
Abell 3667 where the cold fronts are located. The fiber 
allocations for each of these configurations were done au- 
tomatically using the AAOmega specific CONFIGURE 
program 5 . This program uses a "simulated annealing" 
algorithm (jMiszalski et all l2006f ) to ensure the optimal 
allocation of fibers in terms of maximizing the num- 
ber of objects observed. In doing so, it uses the rank- 



ings/weights assigned to targets (§ 12 - 1|) to determine 
what priority they have in being allocated a fiber. 

The use of the simulated annealing algorithm required 
the input catalog be reduced to no more than ~ 800 ob- 
jects for fiber allocation. To achieve this, the catalog was 
first split into the 1452 higher priority objects, on or blue- 
ward of the of the red sequence, and the 711 lower prior- 
ity ones, redward of the red sequence (§ 12. ip . The higher 
priority group were further subdivided into 830 objects 
with r F < 18 and the 622 with 18 < r F < 19. Two fiber 
configurations were required to observe the high prior- 
ity rp < 18 catalog. The small number of objects in the 
r F < 18 catalog that were not observed in these two con- 
figurations were added to the 18 < r F < 19 high priority 
catalog. Finally we constructed a catalog made up of 
those objects not already allocated, plus previously allo- 
cated objects with unreliable redshifts (see below) and a 
small number of low priority objects. Hence in total, six 
configurations were observed, the majority of which tar- 
geted objects lying on or blue-ward of the Abell 3667's 
red sequence. Observing conditions were generally cloud- 
free, but the seeing was rather poor, being in the range 
2.4-4.0 arsec (FWHM). However, given the extended na- 
ture of our galaxy targets due to their relatively low red- 
shifts, mediocre seeing did not have any serious impact 
on our program. The dates, exposure times, magnitudes 
and conditions for the observations are summarized in 
Table [TJ 

For sky subtraction, > 40 fibers per configuration were 
allocated to blank sky regions, selected using the "Gen- 
erate sky grid' command in the CONFIGURE software 
and visually inspected to ensure they were free of objects. 
Approximately 35 known "fringing" fibers were excluded 
from the configuration process. Prior to each fiber con- 
figuration being observed, tungsten and arc lamp expo- 
sures were taken for flat-field and wavelength calibra- 
tion, respectively. The data were fully reduced at the 
telescope using the AAO 2dfdr 6 pipeline reduction sys- 
tem. This system extracts and produces fully flat-fielded, 
wavelength-calibrated and sky-subtracted (but not flux- 
calibrated) spectra for each exposure frame, co-adds the 
spectra in each exposure set, and then splices together 
the blue and red-arm spectra. 

Redshifts were identified and measured from the re- 
duced spectra using the runz code written by Will Suther- 
land for t he Two Degree Fie ld Galaxy Redshift Survey 
(2dFGRS: IC"olless et al.l l2Q 0JJ). This program utili zes the 
cross-correlation method of iTonrv fc Davisl (|1979h . Each 
spectrum was visually inspected and given a redshift 
quality classification, Q, ranging from 1-6, as described 
in Paper I. A Q value of 3, 4, or 5 indicates a reliable red- 
shift was obtained (with a larger value indicating a higher 
quality spectrum) , whereas Q values less than 3 indicate 
either an unreliable redshift or no redshift was obtained. 
A value of Q = 6 was used in cases where the redshift in- 
dicated the object to be a star. We obtained 1714 spectra 
during the run, of which 1570 yielded reliable (Q > 2) 
redshift measurements containing 910 extragalactic ob- 
jects and 660 stellar objects (Q = 6), indicating the un- 
reliability of the SSS star/galaxy classifications. 

2.3. Redshift Measurement Uncertainties 



5 See: http://www.aao.gov.au/AAO/2df/aaomega/aaomega_software.htmtee: http://www.aao.gov.au/AAO/2df/aaomega/aaomega.html 
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TABLE 1 

Summary of the AAT/AAOmega observations for Abell 3667. 



Date 


Object Priority 


Magnitude 


Frames 


Seeing 


2007 Jul 14 


High 


r F < 18.0 


4 X 1800s 


2.8-2.9" 


2007 Jul 16 


High 


r F < 18.0 


3 x 1800s 


2.8-3.5" 


2007 Jul 16 


High 


r F < 19.0 


4 x 1800s 


2.4-3.0" 


2007 Jul 16,18 


High 


r F < 19.0 


2 x 1800 + 3 x 1200s 


3.5" 


2007 Jul 17,18 


High+Low 


r F < 19.0 


3 x 900 + 3 x 1800s 


2.8-4" 


2007 Jul 18 


High+Low 


r F < 19.0 


3 x 1200s 


3.5" 



Of the objects observed twice, 50 had two reliable 
(Q > 2) redshift measurements and so in these cases 
we were able to compare redshifts to assess the robust- 
ness of our measurements and determine their uncer- 
tainty. After excluding one obvious outlier where the 
redshift difference was Acz = 12732km/s, we found 
Acz = 50 ±20 kms -1 for the difference between the first 
and second groups of measurements. We also have 117 
high qualit y measurements in common w ith the redshift 
catalog of Ijohnston-Hollitt et al.1 (|2008l ) and we found 
Acz = — 18 ± 9 kms -1 , again after excluding one outlier 
with Acz = 16134 km s -1 . These redshift comparisons 
are encapsulated in Figure [2j noting that we only con- 
sider objects in the redshift range 0.045 - 0.065 (i.e. that 
appropriate to Abell 3667) . The upper panel of Figure [2] 
shows that the redshift comparisons scatter about a one- 
to-one relationship (solid line with slope unity), indicat- 
ing there are no rcdshift-dependent systematic effects. 
The redshift differences are shown in the lower panel. 
We derive an RMS redshift difference of 151 km s -1 for 
the AAOmega double observations, which implies a red- 
shift uncertainty of 107kms _1 for the sample, compa- 
rable to the RMS of the individual redshift uncertain- 
ties which is 112 km s . An external check is given by 
the RMS of the differences b etween our AAOmega and 
Ijohnston-Hollitt et al.l (|2008h 's redshifts, which gives an 
uncertainty of 64kms _1 , much lower than the quadra- 
ture sum of 148 km s " 1 for the uncertainty abov e and 
the uncertainty of the (|Johnston-Hollitt et al.| [2008) sam- 
ple (~ 100 km s -1 ). We take the precision of our red- 
shift measurements to be the more conservative value of 
107kms" 1 . 

In cases where a galaxy had multiple redshift mea- 
surements, the following policy was adopted in deter- 
mining which value would be used for our analysis. If 
both the measurements were made with AAOmega, we 
adopted the redshift with the higher Q value. If the 
Q values were equal, then we took the redshift with 
the smallest measurement uncertainty (as provided by 
the runz pr ogram). If the second r edshif t measurement 
came from Ijohnsto n-Hollitt et al.l |2008)'s catalog, we 
used our AAOmega redshift (for consistency). Only 
25 extra redshifts wer e adde d to our catalog from the 
Ijohnston-Hollitt et al.l (|2008l ) catalog, taking our total 
number of redshifts to 1595. 

2.4. Redshift Completeness 

The completeness of our redshift measurements was 
determined by calculating the fraction of galaxies in the 
parent photometric catalog for which reliable redshifts 
were obtained. In the left panel of Figure [3] we plot the 
completeness as a function of cluster-centric radius for 
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Fig. 2. — Top panel Shows a comparison of redshift measure- 
ments for multip ly observed obje c ts wh ere the green points show 
AAOmega- Johnston-Hollitt et al. (2008) redshift comparison and 
red points show the comparison of the AAOmega spectra with dou- 
ble observations. The solid black line shows the one-to-one relation- 
ship. The lower panel shows the redshift difference as a function 
of redshift where the color coding is as for the upper panel. 

the magnitude ranges rp < 16, 16 < rp < 17, 17 < 
rp < 18 and 18 < rp < 19. It shows that we have 
achieved excellent completeness within a cluster-centric 
radius of 2.5 Mpc. We also plot the completeness as a 
function of rp magnitude in right-hand panel of Figure 
[3] which demonstrates a high level of completeness across 
all magnitude bins. The cluster is very well sampled 
at magnitudes brighter than rp — 19 for all radii less 
than 2.5 Mpc. I f we assume the cluster M* F = -20.84 
(|Eke et al.ll2004h and neglect K-correction and galactic 
extinction terms, we probe ~ 3 mags down the luminosity 
function. 

3. CLUSTER MEMBER SELECTION 

The major driver of this study is the detection of 
substructure using a sample of spectroscopically con- 
firmed cluster members. To achieve this, it is imper- 
ative that the cluster member catalog be free of line 
of sight interlopers which may produce false substruc- 
ture. Thus, a robust selection procedure is required, 
and we have proceeded in much the same manner as 
described in Paper I. Briefly, initial rejection of fore- 
ground and background objects was performed by sort- 
ing objects in redshift space, determining the velocity 
gap between neighboring galaxies (here velocity is de- 
fi ned as cz) and identifyi ng the cluster (via the method 
of iDe Propris et ai1l2002f) as the peak in redshift space 
that is separated from foreground and background galax- 
ies by a gap greater than some specified value. Here we 
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Fig. 3. — Left: Spectroscopic completeness as a function of cluster-centric radius, for different magnitude ranges: rp < 16 {solid line), 
16 < rp < 17 (dashed line), 17 < rp < 18 (dot-dashed line), and 18 < rp < 19 (dotted line). Right: Spectroscopic completeness as a 
function of rp magnitude. 





Fig. 4. — Redshift distribution for all robust measurements in 
the Abell 3667 field. Abell 3667 clearly stands out against the pro- 
jected background and foreground objects. The arrows indicate the 
separation of the cluster based on gaps in the velocity distribution 
of cz = 1100 km s — 1 (see text). The bin at z = contains some 
660 stellar objects. 

use 110 kms" 1 , which is the velocit y dispersion mea- 
sured bv lJohnston-Hollitt et alj (|2008D . The redshift dis- 
tribution of galaxies within the Abell 3667 field is shown 
in Figure [4] where Abell 3667 stands out clearly from 
foreground and background objects. 

The membership was further refined using the "shifting 
gapper" method of lFadda et al.l (|1996f) which is described 
in detail in Paper I. Briefly, the data are binned radially 
such that each bin contains 50 or more objects and inter- 
lopers are rejected in each radial bin based on gaps in the 
peculiar velocity distribution. Figure [5] shows the results 
of applying the shifting gapper method, where it is clear 
that a number of foreground galaxies are rejected with 
this refined procedure, while a number of background 
galaxies very close in redshift space at radii greater than 
2.5 Mpc are also rejected. 

Our final cluster sample contains 550 members within 
a cluster-centric radius of ~ 3.5 Mpc, more than doubling 
the number of spectroscopically confirmed cluster mem- 
bers previously known. Using this expanded sample, we 
apply the robust biweight scale and locations estimators 



Fig. 5. — Refinement of cluster membership using the shifting 
gapper method (see text for details). The black crosses are cluster 
members, red asterisks indicate objects that are rejected as inter- 
lopers and the blue circles indicate objects rejected prior to the 
shift gapper refinement. 

described in iBeers et"aH Jl990) to determine a cluster 
velocity dispersion, a v = 1056 ± 38 kms -1 , and mean a 
redshift, z cos — 0.0553 ± 0.0002 where the quoted errors 
are ler values and were determined using the jackknifc 
re-sampling technique. Presented in Figure [5] are the 
integral and differential profiles of the location (^t(R)) 
and scale (a(R)) estimators. The differential profiles are 
determined using the radial bins from the shifting gap- 
per analysis. The integral profiles are obtained through 
sorting the members by cluster-centric radius and mea- 
suring fi(R) and cr(i?) within the projected radius of each 
galaxy, beginning with the tenth galaxy from the center. 
The integrated <j(R) becomes relatively constant outside 
R = 2 Mpc, indicating that the cluster velocity disper- 
sion measured at this radius is free from the effects of 
velocity anisotropy, and is representative of the cluster's 
gravitational potential. The differential cr(R) is a declin- 
ing function of radius, and is significantly lower in the 
outskirts. The differential fi(R) measurements are noisy, 
but are generally consistent, within the errors, with the 
value measured for the whole cluster. 

4. SUBSTRUCTURE DETECTION 
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Fig. 6. — Top panel shows the integral (points and solid lcr con- 
fidence limits) and differential (circles with lcr error bars) o~(R) 
measurements. Also plotted are the values measured considering 
all members (solid line with dashed ltr limits), and the value ob- 
tained by assuming the same specific energy in the gas and galaxies 
and using the gas temperature to determine the expected veloc- 
ity dispersion (dot-dashed line with dotted 90% confidence limits). 
Bottom panel shows the corresponding ^i(R) measurements. 

In this section, we search for substructure using the 
velocity distribution, the distribution of galaxies on the 
sky and a combination of spatial plus velocity informa- 
tion. These various statistical methods have differing 
sensitivities to t he variety of substructure that can oc- 
cur i n clusters (|Pinknev et al. 1996; iGirardi fc Bivianol 
2002). Using a range of tests improves our chances of 
detecting any substructure that is present. 

4.1. Quantifying the Gaussianity of the Velocity 
Distribution 

Dynamically relaxed clusters are expected to have line- 
of-sight velocity distributions which are very close to 
Gaussian. Departures from Gaussianity in cluster veloc- 
ity distributions can signify a number of different physi- 
cal phenomena, eg., highly radial orbits, interloper con- 
tamination, circular orbits and, most commonly, sub- 
structure within the cluster. It must be noted, however, 
that non-detection of departures from Ga ussianity is not 
neces sarily evidence of a relaxed system (jPinknev et al.l 
119961 See also Paper I). 

Our initial test for departures from Gaussianity in- 
volves the use of th e Kolmogorov-Smirnov (K-S) statistic 
(|Press et al.lll992l) which measures the largest absolute 
departure between the cumulative probability distribu- 
tion of a data sample and a theoretical distribution de- 
scribing the data (here, a Gaussian) and determines the 
significance of the difference by calculating the likelihood 
that the data are drawn from a parent distribution which 
is Gaussian, taking into account the sample size. The K- 
S test reveals no evidence for non-Gaussianity at the 50% 
level. 

The disadvantages of the K-S test are its lack of sen- 
sitivity in the tails of a distribution (it is most sensitive 
in the vicinity of the median), its spurious results in the 
presence of interlopers and also its failure to quantify the 
way in which a distribution may differ from Gaussianity, 
which is clearly important in determining the physical 
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Fig. 7. — Top panel: peculiar velocity histogram with bins of 
width 200 km s -1 . The solid curve is a Gaussian with a mean of 
zero and standard deviation of 1056 km s — 1 and the dashed curve 
shows the Gauss-Hcrmitc reconstruction of the distribution. Bot- 
tom panel: a strip density plot where the arrows indicate the po- 
sitions of weighted gaps with normalized values greater than 2.5 
(sec text). 



driver of the departure from Gaussianity. Thus, we used 
the metho d first applied to galaxy cluster velocity distri- 
butions bv lZabludoff et all <|1993f) where the distribution 
is approximated by a series of three Gauss-Hermite func- 
tions: the coefficient of the zeroth order term, ho, mul- 
tiplies a best-fitting Gaussian, the third order term, co- 
efficient /13, describes asymmetric deviations from Gaus- 
sianity (similar to skewness) and the fourth order term, 
coefficient /14, describes symmetric deviations (similar 
to kurtosis). It has the advantage of being insensitive 
to outliers in the tails of the distribution (compared to 
the Gaussian skewness and kurtosis measurements). For 
Abell 3667, the third and fourth order terms had values 
of h 3 = 0.009 and /i 4 = 0.047. 

To determine the significance of the results, we gener- 
ated 10,000 Monte-Carlo realizations of Gaussians with 
N = 550 and mean and dispersion equal to the best 
fitting value derived for our data, and determined the 
number of times our observed values occurred in the sim- 
ulations. Values of |ft.3| > 0.009 occur in 77% of the re- 
alizations, so we conclude the observed ft. 3 term is not 
significantly different from zero. Values of \hn\ > 0.047 
occur only 10% of the time, thus we conclude the ob- 
served term is mildly significant. The positive value 
of the hi term means that the velocity distribution has 
slightly longer tails and is more peaked than the best fit- 
ting Gaussian, suggestive of either radial orbits or sub- 
structure. In Figure [7] we plot the binned velocity dis- 
tribution, with both a Gaussian (with a mean of zero 
and a standard deviation of 1056 km s -1 ) and the Gauss- 
Hermite reconstruction curves overlaid. 

We searched for possible segr egation in the velo city dis- 
tribution using weighted gaps (jBeers et al.lll990f) . where 
the velocities are sorted in increasing order and the ith 
velocity gap is given by gi = Vi+i—Vi. The weight for the 
ith gap is Wi — i(N — i) and the weighted gap is defined 
as y/WiQi- The weighted gaps are normalized through di- 
viding by the mid-mean (MM) of the ordered weighted 
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Fig. 8. — The left-hand most panel shows the spatial distribution for the entire cluster member population with smoothed contours (see 
text) drawn at 10galsMpc~ 2 intervals. The positive and negative velocity maps are shown in the middle and right-hand most panels, 
respectively. The contour levels are drawn at intervals of 5 gals Mpc -2 . In each panel, black crosses are drawn at the positions of the two 
dominant cluster galaxies, and the black arc in the south-east quadrant marks the position of the cold front. 

gap distribution given by 



MM 
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where N = 550. Normalized weighted gaps which have 
values larger than 2.5 occur 1.4% of the time when se- 
lecting gaps at random from a Gaussian distribution 
(|Wainer fe Schachtlll978l h thus 8 such gaps are expected 
for Abell 3667. In fact we find 4 such gaps, shown as 
arrows in Figure [7] where the positions of each peculiar 
velocity in the sample is plotted in a strip density plot. 
We conclude that the weighted gap analysis does not pro- 
vide significant evidence for segregations in the velocity 
distribution of Abell 3667. 

4.2. 2D galaxy Distribution 

While the tests for Gaussianity are useful for detecting 
mergers occurring along our line of sight, where the veloc- 
ity distribution is significantly perturbed, they can be in- 
effective in the case of a merger occurring near the plane 
of the sky (jPinkney et al.lfT996l and Paper I). In these 
cases, projected galaxy density maps are useful tools for 
searching for evidence of substructure, but are contami- 
nated by background or foreground interlopers projected 
along the line of sight when redshift information is not 
available. 

Given the presence of significant substructure ap- 
parent in the pr ojecte d galaxy d e nsity maps of 
iProust et all (|1988fl and ISodre et all (|1992D . and the 
non-detection in of significant substructure in the 
maps of spectroscopica l ly co nfirmed cluster members in 
I Johnston- rlollitt et all (|2008h . we use our much larger 
sample to produce a galaxy density map to search for 
signs of substructure. We plot the distribution of all 
galaxies in the sample, along with the associated sur- 
face density contours in the left-hand panel in Figure [51 
The contours were generated by smoothing the distribu- 
tion with a variable-width Gaussian, where a varies from 
140 kpc in the cen ter to 500 kpc i n the outski rts. We find 
similar results to IProust et al.1 (|1988j ) and ISodre etlrtl 



(1992), finding Abell 3667's galaxy surface density dis- 
tribution is clearly bimodal, with a high surface density 
core at the center and also a high-density region offset 
by ~ 1 Mpc to the north-west which is coincident with 
the second dominant cluster galaxy. There are two other 
clear peaks in the surface density, one ~ 1.2 Mpc to the 
south-east and one ~ 1.1 Mpc just west of due north. 

As a first order test for correlation between multi- 
peaked surface density distribution and the velocity dis- 
tribution, we split the catalog into positive and negative 
peculiar velocities and plot the spatial distributions. The 
two segregations are shown in the middle and right-hand 
panels of Figure [SJ along with the associated galaxy sur- 
face density contours, generated using the same Gaussian 
smoothing described above. The most notable features 
upon comparing the positive and negative maps are; the 
south-east and north-east peaks, which are clearly more 
significant in the negative peculiar velocity map and the 
north-west peak appears in both the positive and nega- 
tive maps, although it is more significant in the positive 
velocity map and slightly offset in the negative velocity 
map. 

4.3. 3-D tests for sub-structure 

In the previous section, we searched for spatial/velocity 
correlations by crudely splitting the catalog into posi- 
tive and negative velocity samples. In this section, we 
combine both velocity and spatial information to test 
for substructure in a more refined manner. In general 
it has been found that tests simultaneously employing 
both spatial and velocity inf ormation are most su ccessful 
at detecting substructures (Pin knev et al.|[l996l ). Here, 
we incorporate the k-statistic, hereafter k s , which was 
fir st used in the analysis of substructure within Coma 
by iColless fe Durml (|1996f ). The k s tests for departures 
between local and global kinematics by determining the 
probability that the local velocity distribution differs 
from the global one. In brief, the evaluation of k s is 
performed as follows: The nearest n = ^/N neighbors in 
projected radius from each galaxy in the sample are se- 
lected, where N is the number of cluster members. The 
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velocity distribution of the n neighbors is compared to 
the velocity distribution of the remaining N — n galaxies 
using the K-S D-statistic and the likelihood of finding a 
D-statistic larger than the observed D-statistic is deter- 
mined. Hence, low likelihoods give an indication that the 
local and global velocity distributions are significantly 
different. The negative log of the likelihoods are summed 
to give an overall measure of the substructure present 
within the cluster. The significance of the substructure 
is determined by repeating the above analysis for 10,000 
Monte Carlo realizations, where any correlation between 
the spatial and velocity distributions due to substructure 
within the cluster is erased by randomly shuffling the pe- 
culiar velocities in the sample, whilst maintaining spatial 
information. For Abell 3667, we measured an overall k s 
of 422. In contrast, the k s distribution produced in our 
Monte Carlo realizations is well modeled by a log-normal 
distribution with a mean of )i(hiK s )= 5.44 and standard 
deviation cr(ln k s )=0.15, indicating our result is signifi- 
cant at the 3.96a level, and a value as large as 422 was 
n ot obtained in a ny re alization. 

iPinknev et al.1 (|1996l ) showed that strong radial gradi- 
ents in the cluster velocity dispersion profiles can pro- 
duce high false positive detection rates for 3D substruc- 
ture tests. With this in mind, and given that Fig- 
ure [6] shows a significant gradient in the a(R) profiles 
for Abell 3667, we re-measured k s within the 2Mpc re- 
gion shown in Figure [3J The observed k s in the smaller 
region is 312, and a value of this order is obtained in 
only 12 of our 10,000 Monte Carlo realizations which 
produce a mean of /i(ln k s )= 5.27 and standard deviation 
cr(ln k s )=0.17, indicating the observed k s is significant 
at the 3.13cr level, again indicating the existence of sub- 
structure within Abell 3667 at high significance. 

The most effective way to illustrate the substructure 
within the cluster is to present 'bubble' plots, with cir- 
cles whose radii proportional to the negative log of the 
measured K-S likelihood are plotted at each cluster mem- 
ber position. This is shown in Figure [9l where regions of 
clustered large circles indicate the existence of substruc- 
ture. Circles in this plot are color-coded, red for galaxies 
with positive v pec , and blue for galaxies with negative 
v pec . Also, galaxies having K-S likelihoods which occur 
in only 5% of the simulated cases are emboldened. The 
circles are plotted on top of the galaxy surface density 
contours, generated as described in ^4.21 

There are a number of regions where the presence of 
groupings of significantly large 'bubbles' indicate local- 
ized velocity structure, confirming our measure of sig- 
nificant substructure being present within Abell 3667. 
These localized regions of velocity structure are gener- 
ally coincident with surface density enhancements, and, 
most significantly in terms of the cold front, we observe 
a galaxy density enhancement 700 kpc to the south-east 
of the front which is coincident with a clustering of large 
bubbles, indicating the presence of a subgroup. Interest- 
ingly, the significant enhancement in galaxy density to 
the north-west, coincident with the 2nd ranked cluster 
galaxy, shows little to no deviation in its local veloc- 
ity distribution, despite bein g coincident with a peak in 
the projected mass maps of Uoffre et al.l (|2OOO0 . There 
is, however, localized velocity structure surrounding the 
substructure which may be associated. 

It is possible that the non-detection of localized ve- 
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Fig. 9. — A 'bubble' plot of Abell 3667, where the circles rep- 
resent the size of the k-statistic, which provides a measure of the 
departure of the local velocity distribution from the global one. 
The circles are colored blue or red, depending on whether v pec is 
negative or positive, respectively. Circles plotted with a thicker 
line indicate a significant departure whereby the k-statistic value 
occurs only in 5% of 10,000 Monte Carlo realizations. Green cir- 
cles indicate regions of interest which are used as inputs to the 
KMM algorithm, and the black dot-dash circle shows the region 
within which the KMM analy sis w as done. The black contours 
are generated as described in § 14.21 and are plotted at intervals of 
10 gals Mpc -2 and black crosses are drawn at the positions of the 
two dominant cluster galaxies, and the black arc in the south-east 
quadrant marks the position of the cold front. 



locity structure around the region of enhanced galaxy 
density to the north-west is due to the difference in the 
scale probed by the test compared to the physcial scale 
of the substructure. Since the local velocity distribu- 
tion is defined using a fixed number of near neighbors, 
rather than the near neighbors within some fixed pro- 
jected radius, it is difficult to determine exactly which 
scales are being probed at each galaxy position in Fig- 
ure El All that is known is that where the galaxy density 
is high, smaller scales are probed and, conversely, where 
the galaxy density is low, larger scales are probed. To 
remedy this, and also to address the fact that substruc- 
tures can have a range of different physical scales, we use 
a multi-scale analysis. This involves remeasuring k s in a 
similar fashion to that described above, but rather than 
using a fixed number of n nearest neighbors to define the 
local velocity distributions, we use the near neighbors 
within a fixed projected radius, Rnn, of the galaxy of in- 
terest. The analysis is repeated for near neighbors within 
Rnn — 200, 400, 800 and 1600 kpc and the resulting bub- 
ble plots are shown in Figure[l0] For these radii, we mea- 
sure k s values of 342 (4.15cr), 405 (3.71er), 762 (4.62a) and 
1114 (4.04er), respectively, where the bracketed values are 
the significances, measured as above from 10,000 Monte 
Carlo realizations. Across the four scales presented in 
Figure |T0l all of the substructures seen in Figure [9] are 
reproduced and, most significantly, on scales 800 and 
1600 kpc the north-west substructure appears as a region 
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Fig. 10. — Multi-scale bubble plots where the circle size and color is defined as in Figure [9] but rather than using the n = \HV nearest 
neighbors to define the local velocity distribution, we use the near neighbors within a fixed projected aperture of i?jvjv=200, 400, 800 and 
1600 kpc shown in the (upper left), (upper right), (lower left) and (lower right) panels, respectively. The bold black circle in the upper left 
of each panel has radius -Rjvjv- 



containing significant local velocity substructure. 

Our concern regarding the <r(R) gradient affecting the 
k s measurement seems warranted, given the detection of 
significant velocity structure > 2 Mpc to the south-east 
in Figure[5]and the lower panels of Figure fTUl This struc- 
ture is not offset significantly in peculiar velocity, how- 
ever its dispersion is much lower than the total cluster 
value and therefore it is unlikely to be a bound infalling 
substructure. Rather it is likely to be an artifact of the 
declining a(R) profile. 

Having identified significant substructure within 
Abell 3667, we now use the KMM (Ka ye's Mixture 
Modeling) algorithm ([Ashman et al.lfl994( ) to partition 
the different substructures. The algorithm fits a user- 
specified number of Gaussians to a data set and deter- 
mines the improvement of the fit over that of a single 
Gaussian. Here we utilize the combination of spatial and 



velocity information and fit 3-D Gaussians. Inspection 
of Figure [5] reveals four regions of interest (outlined by 
green circles) where there is localized velocity structure 
and/or there is structure in the galaxy surface density. 
These regions serve as an initial estimate of the number 
of partitions into which the KMM algorithm should di- 
vide the cluster, and we restrict our analysis to the region 
containing these structures defined by the black dot-dash 
region in Figure [5] Within the regions of interest, ob- 
vious velocity outlyers are rejected, and the mean and 
standard deviation are measured for the spatial distribu- 
tions, while the m edian and median absolute deviation 
([Beers et al.l Il990h are measured for the velocity distri- 
bution. The measurements serve as initial inputs to the 
algorithm and, along with the outputs from the KMM 
algorithm, are shown in Table [2J where the final v and 
a v values are determined using the median and gapper 
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estimators in cases where Nq < 15, and the biweight 
location and scale estimators in cases where Ng > 15 
members, respectively. 

In Figure 1111 we show the spatial distribution of the 
partitions, along with the associated velocity histograms 
for each partition. The cluster is essentially divided into 
north-west and south-east partitions, with the south-east 
one associated with the dominant cluster galaxy, and the 
north-west one associated with the second ranked galaxy 
coincident with the second major overdensity in the 
galaxy surface brightness maps. The velocity dispersions 
of these components are very similar (1039 ± 66kms _1 
for KMM2 and 1073 ± 64kms- 1 for KMM5), but their 
velocities are slightly offset, with KMM2 having v = 
422 ± 82kms~ 1 , similar to the peculiar velocity of the 
of the second ranked galaxy which is at 432kms _1 , and 
KMM5 has v = -103 ± 71 kms" 1 . There are two other 
minor structures (KMM1 and KMM3) which appear to 
be small bound groups, and a low velocity dispersion 
structure to the south-east (KMM4) , which is associated 
with an overdensity in the galaxy surface density maps. 

5. DISCUSSION AND CONCLUSIONS 

Abell 3667 shows three signatures of a recent merger. 
There is significant substructure in both the spatial and 
localized velocity distributions of the galaxies. It has a 
pronounced c old front and other structure in its X-ray 
emitting gas (iKnopp et al.lll996UMarkevitch et alJH99a 
IVikhlinin et al. 200 lbl) and it contains two radio relics 
(jRottgering et al.l Il997). Here we consider how these 
structures were created. 

5.1. Merger Scenario 

The multiwavelength data at hand may be reconciled 
with a simple two-cluster merger. The structure to the 
north-west is clearly significant in the 2D galaxy density 
maps and is coincident with the second dominant clus- 
ter galaxy. The structure is also coincident with a sig- 
nificant mass conc entration in the weak-lensing maps of 
iJoffre et al.l ((2000). There is only a small offset in pecu- 
liar velocity of ~ 422 km s -1 , and very little difference in 
the velocity dispersion between the north-west substruc- 
ture and the main cluster. This suggests a merger be- 
tween subclusters of similar masses taking place roughly 
in the plane of the sky. 

The direction of motion of the north-west substruc- 
ture cannot be determined from the optical data alone, 
however it can be constrained by considering the radio 
and X-ray data in the literature. Regarding the radio 
data, the pertine nt observation is the nor th-west radio 
relic described in IRottgering et al.l (|1997h where infer- 
ence can be made regarding the direction of motion by 
consideration of the s pectral index and brigh tness dis- 
tribution of the relic. IRottgering et al.l (1997) find that 
the relic brightness steepens at the north-west edge and 
is more diffuse to the south-east, while the spectral in- 
dex is flat at the north-west edge and steepe n s tow ards 
the south-east. The model of iRoettiger et al.l ()1999l ) ex- 
plains these features as due to a shock front moving from 
the south-east to north-west, where the steepening of the 
brightness and the flatness of the spectral index at the 
north-west edge occur because this is where the shock 
acceleration of electrons is currently occurring. Since 
the shock moves from south-east to north-west, the elec- 
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Fig. 11. — Upper panel: The spatial distribution of the KMM 
partitions. The ellipses are the 2<r boundaries for the Gaussians 
fitted to the spatial distributions. Galaxy contours with the same 
spacings as in Figure [9] are shown in pink, and the color key for the 
positions is given at upper left. Black crosses are drawn at the posi- 
tions of the two dominant cluster galaxies, and the black arc in the 
south-east quadrant marks the position of the cold front. Bottom 
panel: the velocity histograms corresponding to the KMM parti- 
tions, color-coded to match the upper panel. At the bottom right 
we show the velocity distribution for all galaxies included in the 
KMM analysis, along with the break-down of the Gaussian compo- 
nents making up the total velocity distribution. We also ovcrplot 
the best fitting single Gaussian {pink curve) for comparison. 

trons to the south-east no longer have a source of en- 
ergy for acceleration, and the spectral index steepening 
is caused by the aging of the particles as the shock moves 
off to the north-west. It should be noted, ho wever, that 
recent observations (jJohnston-Hollittl l2003h have cast 
doubt on the existence of a flat spectral index at the 
north-west edge of the relic and show only a moderate 
spectral gradient across the source. Even without these 
arguments, the curvature of the radio relic is consistent 
with a shock propaga ting to the northwe st. The ROSAT 
X-ray observations of lKnopp et al.l (jT996) find significant 
evidence for extended X-ray emission associated with the 
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TABLE 2 

Input and output parameters for the KMM algorithm. The final column gives the estimates for the 

CORRECT ALLOCATION RATE FOR EACH KMM PARTITION. 





Initial Inputs 






KMM Outputs 






Group 










{(T x , 0-y,(T v ) 


N S al 


Rate (%) 


KMM1 


(130, 1145, -1388) 


(112, 66, 353 ) 


10 


(153, 1169, -1483) 


(83, 58, 395) 


9 


100 


KMM2 


(905, 747, 408) 


(368, 370, 829) 


110 


(879, 824, 422) 


(530, 560, 1039) 


164 


89 


KMM3 


(528, -373, 1375) 


(91, 89, 704) 


9 


(563, -415, 1457) 


(64, 44, 426) 


7 


100 


KMM4 


(-542, -798, -523) 


(292, 238, 489) 


36 


(-761, -659, -638) 


(412, 326, 209) 


27 


73 


KMM5 


(0, 0, 55) 


(827, 853, 1020) 


284 


( -144, -173, -103) 


(676, 739, 1073) 


242 


87 


Note. - 


— The units of aJ, y , a x 


and (7 y are kpc, 


and the units of v and a v arc 


kms -1 . For T; and 


(T v , the median and 



gappcr estimators are used when Ng < 15 and where Nq > 15, the biweight location and scale estimators arc used. 



north-west subcluster. In Figure [T2l we present the Dig- 
ital Sky Survey r-band image of Abell 3667 with X-ray 
contours from Chandra and ROSAT overlaid. The con- 
tours clearly extend to encompass the north-west sub- 
structure. The fact that the subcluster hosts an X-ray 
halo, together with the arguments based on the radio 
data, suggest that the north-west subcluster has signifi- 
cant mass and is moving from south-east to north-west 
with the majority of its motion in the plane of the sky. 

In Figure [12] the cold front is seen as a compression 
of the X-ray contours t o the so uth-east of the cluster 
center. iVikhlinin et al.1 (|2001aft argued that the cold 
front is moving in a south-easterly direction at approx- 
imately the sound speed of the ICM. Based on argu- 
ments about the constant shap e and the stability of the 
front, |Vikhlinin & Ma rkevitchl (|2002l ) predict the exis- 
tence of a dark halo associated with the front, with its 
centroid ~ 90 kpc from the leading edge of the front. 
We find no evidence for any structure in the galaxy dis- 
tribution in this region, and n ote that the weak lens- 
ing maps of Uoffre et al.1 (|2000f) do not show any signif- 
icant mass concentrations here either. In the two-body 
merger scenario, the cold front is produced when core 
gas is 'sloshed' out of the central potential well after 
the core passage of KMM2, a nd an excellent example 
of this is seen in iPoole et all (pool 's Fi gure 18 which 
shows entropy maps from the 3:1 mass ratio, offset im- 
pact parameter_shnulations. Interestingly, the structure 
that iMazzotta et al.l (|2002h interpret as a hydrodynamic 
instability is seen in these simulations and interpreted as 
a 'plume' of cool gas ejected from the main cluster's cool 
core due to the subcluster's pericentric passage. These 
simulations also produce two shock fronts on opposite 
sides of the main cluster moving in opposite directions, 
which would account for the observed radio relics in 
Abell 3667. This scenario is essen tially the same as that 
propose d inlRoettiger et all (|1999D . although the simula- 
tions of IPoole et alJ (|2006) are at higher resolution and 
thus provide a better view of the cold front. In this sce- 
nario, KMM4 can be interpreted in two ways; In the first, 
KMM4 is a group either infalling from the foreground, 
or is a background structure that is not physically asso- 
ciated with the cluster, nor with any of the X-ray and 
radio features or the merger. The second interpretation 
is that KMM4 was once associated with KMM2, and was 
tidally stripped from the main KMM2 structure during 
its approach to the main cluster core and survives as a 
dynamically distinct substructure within the cluster. 

An alternative scenario for the features seen in 



Abell 3667 involves a three body merger, between 
KMM2, KMM4 and KMM5 where the arguments re- 
garding KMM4 and its association with the X-ray and 
radio structures observed in the north-west presented 
above are essentially unchanged, however the cold front 
is the remnant core of KMM4, the centroid of which 
lies ~ 550 kpc to the south-east of the front. A shock 
front driven by KMM4 is responsible for the radio relic 
in the south-east. In this scenario KMM4 has traveled 
from the north-west and passed through the cluster core 
where its ICM was significantly slowed by ram pressure, 
while the collisionless galaxies and dark matter continued 
their passage towards the south-east, similar to that seen 
in the B ullet cluster dMarkevitch et a l. 2002). The mass 
maps of I Joffre et al.1 (|2000l ) indicate a significant mass 
concentration not quite coincident with KMM4, but off- 
set by ~ 300 kpc to the south-west. This lies well within 
the 2a boundary seen in Figure fTTl and the mass is highly 
likely to be associated with KMM4. The velocity disper- 
sion of KMM4 is low and is unlikely to reflect the ini- 
tial mass of the system, since the group must have been 
strongly affected by tidal forces and dynamical friction 
during its passage through the cluster core. The mass of 
the group was probably not comparable to the Abell 3667 
system, since its gravitational potential well was unable 
to stop th e gas from being completely str ipped from the 
group (see lAscasibar fc Markevit ch 2006, for a compari- 
son of the effects of different mass ratios on ram pressure 
stripping of core gas). 

5.2. Summary and Conclusions 

We have presented an optical analysis of the cold front 
cluster Abell 3667 based on a large sample of (550) spec- 
troscopically confirmed cluster members using data ob- 
tained on the AAT AAOmega facility. We find the clus- 
ter has a velocity dispersion of a = 1056 km s -1 and 
redshift z = 0.0553. Although the velocity distribution 
does not differ significantly from a Gaussian, the clus- 
ter shows clear evidence of subclustering in the galaxy 
density plots, and also contains significant local veloc- 
ity substructure according to the k-statistic. Using the 
KMM algorithm, the cluster can be partitioned into 
3 main components - the main cluster (KMM5) with 
242 members, a = 1073kms~ 1 , v pec = — 103kms _1 , 
a north-west subcluster (KMM2) with 164 members, 
a = 1039 km s -1 , V pec — 422 km s -1 , and a south-east 
subgroup (KMM4) with 27 members, a = 209 km s -1 
and Vp ec — — 638 km s -1 . 

We confirm that Abell 3667 is a merging cluster, 
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Fig. 12. — Digital Sky Survey r-band image of Abell 3667. Circles mark the cluster members allocated to groups using the KMM method 
and are color coded to match the key seen in the upper panel of Figure [TT1 The black contours are from archival Chandra (central 5 contour 
levels) and ROSAT (outer 5 contour levels) X-ray images. The cold front is visible as a com pression of th e X-ray contours to the south-east 
of the cluster center. The significant excess of emission associated with KMM2 noted by Knopp et al. (1996) is visible as an elongation 
of the contours toward the north-west, encompassing the KMM2 substructure. Black squares mark the positions of the two domninant 
cluster galaxies. 



and the physical phenomena associated with observa- 
tions at X-ray and radio wavelengths can be associ- 
ated with substructure observed in the galaxy distribu- 
tion, although it is apparent that detailed, combined N- 
body/hydrodynamical simulations are required to deter- 
mine the plausibility of the two scenarios presented. This 
study highlights the need for large, deep spectroscopic 
samples to be used in conjunction with the full multi- 
wavelength observations when determining the dynam- 
ical state of a cluster of galaxies. Future work will in- 
volve using the spectroscopic sample to study the galaxy 
properties within Abell 3667, and to correlate them with 
substructures within the cluster. 
We have now studied in detail two clusters from our 



sample of 'cold front' clusters selected from the Chandra 
archive, finding here, as for Abell 1201 in Paper I, that 
despite difficulties in determining precise merger histo- 
ries, cold fronts can be directly related to cluster merger 
activity and thus are reliable signposts of ongoing merg- 
ers. However, it must be noted that both Abell 1201 
and Abell 3667 harbor signs of merger activity (other 
than cold fronts) in their Chandra images, thus were 
expected to show evidences of merger activity at other 
wavelengths. The acid test for this claim will come from 
detailed observations of cold front clusters with relaxed 
X-ray morphologies. 

6. ACKNOWLEDGMENTS 



Substructure in Abell 3667 



13 



Wc thank David Woods and Gregory Poole for useful 
conversations. We also thank the staff at the Anglo- 
Australian Observatory for their support during the ob- 
servations, and in particular Matthew Colless for offering 
Director's Time in order for this project to be under- 
taken. MSO acknowledges the support of an Australian 
Postgraduate Award. We acknowledge the financial sup- 
port of the Australian Research Council (via its Discov- 



Abell, G. O., Corwin, Jr., H. G., & Olowin, R. P. 1989, ApJS, 70, 
1 

Ascasibar, Y., & Markevitch, M. 2006, ApJ, 650, 102 
Ashman, K. M., Bird, C. M., & Zepf, S. E. 1994, AJ, 108, 2348 
Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 
Bliton, M., Rizza, E., Burns, J. O., Owen, F. N., & Ledlow, M. J. 

1998, MNRAS, 301, 609 
Caldwell, N., & Rose, J. A. 1997, AJ, 113, 492 
Churazov, E., & Inogamov, N. 2004, MNRAS, 350, L52 
Colless, M., Dalton, G., Maddox, S., Sutherland, W., Norberg, P., 
Cole, S., Bland-Hawthorn, J., Bridges, T., Cannon, R., Collins, 
C, Couch, W., Cross, N., Deeley, K., De Propris, R., Driver, 
S. P., Efstathiou, G., Ellis, R. S., Frenk, C. S., Glazebrook, K., 
Jackson, C, Lahav, O., Lewis, I., Lumsden, S., Madgwick, D., 
Peacock, J. A., Peterson, B. A., Price, I., Seaborne, M., & Taylor, 
K. 2001, MNRAS, 328, 1039 
Colless, M., & Dunn, A. M. 1996, ApJ, 458, 435 
De Propris, R., Couch, W. J., Colless, M., Dalton, G. B., Collins, 
C, Baugh, C. M., Bland-Hawthorn, J., Bridges, T., Cannon, R., 
Cole, S., Cross, N., Deeley, K., Driver, S. P., Efstathiou, G., 
Ellis, R. S., Frenk, C. S., Glazebrook, K., Jackson, C, Lahav, 
O., Lewis, I., Lumsden, S., Maddox, S., Madgwick, D., Moody, 
S., Norberg, P., Peacock, J. A., Percival, W., Peterson, B. A., 
Sutherland, W., & Taylor, K. 2002, MNRAS, 329, 87 
Ebeling, H., Voges, W., Bohringer, H., Edge, A. C, Huchra, J. P., 

& Briel, U. G. 1996, MNRAS, 281, 799 
Eke, V. R., Frenk, C. S., Baugh, C. M., Cole, S., Norberg, P., 
Peacock, J. A., Baldry, I. K., Bland-Hawthorn, J., Bridges, T., 
Cannon, R., Colless, M., Collins, C, Couch, W., Dalton, G., dc 
Propris, R., Driver, S. P., Efstathiou, G., Ellis, R. S., Glazebrook, 
K., Jackson, C. A., Lahav, O., Lewis, I., Lumsden, S., Maddox, 
S. J., Madgwick, D., Peterson, B. A., Sutherland, W., & Taylor, 
K. 2004, MNRAS, 355, 769 
Fadda, D., Girardi, M., Giuricin, G., Mardirossian, F., & Mezzetti, 

M. 1996, ApJ, 473, 670 
Girardi, M., & Biviano, A. 2002, in Astrophysics and Space Science 
Library, Vol. 272, Merging Processes in Galaxy Clusters, ed. 
L. Feretti, I. M. Gioia, k, G. Giovannini, 39-77 
Girardi, M., Giuricin, G., Mardirossian, F., Mezzetti, M., & 

Boschin, W. 1998, ApJ, 505, 74 
Hocft, M., & Briiggen, M. 2007, MNRAS, 375, 77 
Joffre, M., Fischer, P., Frieman, J., McKay, T., Mohr, J. J., Nichol, 
R. C, Johnston, D., Sheldon, E., & Bernstein, G. 2000, ApJ, 534, 
L131 

Johnston-Hollitt, M. 2003, Ph.D. Thesis 

Johnston-Hollitt, M., Hunstead, R. W., & Corbett, E. 2008, A&A, 
479, 1 

Knopp, G. P., Henry, J. P., k, Briel, U. G. 1996, ApJ, 472, 125 
Markevitch, M., Gonzalez, A. H., David, L., Vikhlinin, A., Murray, 

S., Forman, W., Jones, C, & Tucker, W. 2002, ApJ, 567, L27 
Markevitch, M., Sarazin, C. L., & Vikhlinin, A. 1999, ApJ, 521, 

526 

Mazzotta, P., Fusco-Femiano, R., & Vikhlinin, A. 2002, ApJ, 569, 
L31 

Miszalski, B., Shortridge, K., Saunders, W., Parker, Q. A., & 
Croom, S. M. 2006, MNRAS, 371, 1537 



ery Project Scheme) throughout the course of this work. 
PEJN was supported by NASA grant NAS8-01130. 

This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. 



Owers, M. S. 2008, Ph.D. Thesis 

Owers, M. S., Nulsen, P. E. J., Couch, W. J., Markevitch, M., & 
Poole, G. B. 2009, ApJ, in press (astro-ph/0810.4650; Paper I) 

Pinkney, J., Roettiger, K., Burns, J. O., & Bird, C. M. 1996, ApJS, 
104, 1 

Poole, G. B., Fardal, M. A., Babul, A., McCarthy, I. G., Quinn, 

T., & Wadsley, J. 2006, MNRAS, 373, 881 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, 

B. P. 1992, Numerical recipes in FORTRAN. The art of scientific 

computing (Cambridge: University Press, — cl992, 2nd cd.) 
Proust, D., Mazure, A., Sodre, L., Capclato, H., & Lund, G. 1988, 

A&AS, 72, 415 

Ramclla, M., Biviano, A., Pisani, A., Varela, J., Bettoni, D., Couch, 
W. J., D'Onofrio, M., Dressier, A., Fasano, G., Kj0rgaard, P., 
Moles, M., Pignatelli, E., & Poggianti, B. M. 2007, A&A, 470, 
39 

Roettiger, K., Burns, J. O., & Stone, J. M. 1999, ApJ, 518, 603 

Rottgering, H. J. A., Wieringa, M. H., Hunstead, R. W., & Ekers, 
R. D. 1997, MNRAS, 290, 577 

Saunders, W., Bridges, T., Gillingham, P., Haynes, R., Smith, 
G. A., Whittard, J. D., Churilov, V., Lankshear, A., Croom, 
S., Jones, D., & Boshuizen, C. 2004, in Presented at the Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference, 
Vol. 5492, Ground-based Instrumentation for Astronomy. Edited 
by Alan F. M. Moorwood and lye Masanori. Proceedings of the 
SPIE, Volume 5492, pp. 389-400 (2004)., ed. A. F. M. Moorwood 
& M. Iye, 389-400 

Sharp, R., Saunders, W., Smith, G., Churilov, V., Correll, D., 
Dawson, J., Farrel, T., Frost, G., Haynes, R., Heald, R., 
Lankshear, A., Mayfield, D., Waller, L., & Whittard, D. 2006, 
in Presented at the Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference, Vol. 6269, Ground-based and 
Airborne Instrumentation for Astronomy. Edited by McLean, 
Ian S.; Iye, Masanori. Proceedings of the SPIE, Volume 6269, 
pp. 62690G (2006). 

Smith, G. A., Saunders, W., Bridges, T., Churilov, V., Lankshear, 
A., Dawson, J., Correll, D., Waller, L., Haynes, R., & 
Frost, G. 2004, in Presented at the Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference, Vol. 
5492, Ground-based Instrumentation for Astronomy. Edited by 
Alan F. M. Moorwood and Iye Masanori. Proceedings of the 
SPIE, Volume 5492, pp. 410-420 (2004)., ed. A. F. M. Moorwood 
& M. Iye, 410-420 

Sodre, L. J., Capelato, H. V., Steiner, J. E., Proust, D., & Mazure, 
A. 1992, MNRAS, 259, 233 

Tonry, J., & Davis, M. 1979, AJ, 84, 1511 

Vikhlinin, A., Markevitch, M., & Murray, S. S. 2001a, ApJ, 551, 
160 

— . 2001b, ApJ, 549, L47 

Vikhlinin, A. A., & Markevitch, M. L. 2002, Astronomy Letters, 
28, 495 

Wainer, H., & Schacht, S. 1978, PsychometriKa, 43, 203 
Zabludoff, A. I., Franx, M., & Geller, M. J. 1993, ApJ, 419, 47 



